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ABSTRACT 

Using near infrared UKIDSS Galactic Plane Survey data, we make extinction mea- 
surements to individual stars along the same line of sight as molecular clouds. Using 
an existing 3D extinction map of the inner Galaxy, that provides line of sight specific 
extinction-distance relationships, we convert the measured extinction of molecular 
clouds to a corresponding distance. These distances are derived independently from 
kinematic methods, typically used to derive distances to molecular clouds, and as such 
they have no near/far ambiguity. The near/far distance ambiguity has been resolved 
for 27 clouds and distances have been derived to 20 clouds. The results are found to 
be in good agreement with kinematic measurements to molecular clouds where the 
ambiguity has already been resolved, using HI self-absorption techniques. 
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1 INTRODUCTION 

Molecular clouds are the birth places of massive stars. 
They form in an environment shielded from Galactic radia- 
tion and so their cold temperatures allow the formation of 
H2. Star formation is triggered through gravitational insta- 
bility inside the cloud, when a cloud fragment's gravitational 
energy exceeds its internal energy. As massive star forming 
regions are retentive of their natal cloud, until it is dispersed 
by the formation of HII regions, if the distance is known to 
the molecular cloud, then the distance is known to the mas- 
sive star forming region. 

The Red MSX Source (RMS) survey (|Urauhart et all 

[2003) has carried out a series of multi-wavelength follow 
up observations, to identify genuine massive young stel- 
lar objects (MYSOs), and ultra compact (UC) HII regions. 
Initially, a large sample of MYSO candidates were colour- 
select ed from the MSX point s ource catalogue (jEgan et al.l 
I2OO3I ) bv lLumsden et al.l (|2002l '). To confirm the presence of 
a massive young star, we need to determine the luminosity 
of each source. Thus, we need an accurate measurement of 
each MYSO candidate's heliocentric distance. 

Every MYSO candidate in the RMS survey has had 
a kinematic distance determined, from ^^CO observations 
lUrguhart et al.l 120081 ) . to the molecular cloud from which 
it was formed. Kinematic distances are derived by applying 
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a Galactic rotational model to measured kinematic veloc- 
ities. Excluding localised velocity perturbations, the kine- 
matic velocity is the projection of the orbital velocity of a 
molecular cloud, about the Galactic centre, along the line of 
sight. Therefore for clouds within the Solar circle, there is 
not a unique solution to the derived distance and such clouds 
will possess a near/far kinematic ambiguity. Urquhart et al 
in prep.) have used a r ecent catalogue of molecular clouds 
Rathborne et al.l [2009|) to identify the molecular clouds 
most likely associated with RMS sources. Using HI self- 
absorption and also absorption towar ds 21 cm continuum 
sources embedded in molecular clouds, iRoman-Duval et al.l 
((20091) resolved the near /far ambiguity to ~90% of all molec- 
ular clouds identified bvlRathborne et al.l (|2009r ). Using data 
from iRoman-Duval et al.l ( 2009f ). Urquhart et al (in prep.) 
have resolved the kinematic ambiguity to 186 RMS sources. 
Kinem atic distance ambigu i ties, although shown to be 
solvable bv lRoman-Duval et al.l ((20091), are not the only po- 
tential pitfall when attempting to derive distances to molec- 
ular clouds using kinematic techniques. The conversion of 
kinematic velocities into kinematic distances requires an ac- 
curate Galactic rotation model, and the assumption that 
the measured kinematic velocity is due entirely from rota- 
tion about the Galactic centre. Kinematic velocities infer 
distances of ~4 kpc to stars found in certain portions of 
the Perseus spiral arm. However there is a discrepancy of 
a fa ctor of two using luminosity distance estimates, ~2.2 
kpc ((Humphrevdl 19781), and by tri gonometric parallax to 4 
masers, 1.95±0.04 ( Xu et al.ll2006l ). Further disagreements 



have been found to 18 masers bv lReid et alj (|2009l ). but par- 
allax measurements are still sparse. Clearly there is much to 
gain from any distance determination method that is inde- 
pendent to Galactic rotation models. 

Molecular clouds are very optically thick. As such, stars 
behind molecular clouds will suffer much more extinction 
relative to stars in front of molecular clouds along the same 
line of sight. Several authors have used extinction mea- 
surements to perform a wide varie ty of tasks. Us i ng cu - 
mulative star counts from 2MASS, iFroebrich et all (120051 ) 
have c onstructed relative extinc tion maps of the Galactic 
Plane. iLombardi fc Alvesl |200j) deredden 2MASS data to 
determine the near infrared colour excess of stars to pro- 
duce extinction maps. They deredden all stars to a sin- 
gle, average intrinsic co lo ur using t h e red dening vector of 
iRieke fc Lebofskvl (|l985l ). ISale et al.1 (|2009l ) use IPHAS Ha 
data to simultaneously determine extinction, intrinsic colour 
and distance estimates of early-A to K4 stars. They use these 
data to map extincti on in three dimensions a cross the north- 
ern Galactic Plane. [Marshall et al.l (|200(i ) also construct 
three dimensional extinction maps of the inner Galaxy using 
2MASS data and the Stellar Population Synthesis Model of 
the Galaxy, developed in Besangon. They used the synthetic 
Besangon data to provide the intrinsic colour and probable 
distances of colour selected giant stars, mapping the inner 
Galaxy in 15'xl5' tiles. 

Using the most reliable United Kingdom Infrared Deep 
Sky Survey (UKIDSS) Galactic Plane Survey (GPS) data, 
this work dereddens individual giant stars along a particular 
reddening track, to the point of intersection on an intrinsic 
giant locus. The dereddening process allows us to measure 
the extinction suffered by stars along the same line of sight as 
molecular clouds. We then determine distances, that are in- 
dependent to Galactic rotation models, to molecular clouds 
using line of sight specific extinction-distance relationships, 
derived from an existing 3D extinction map of the inner 
Galaxy. We acq uire the 3D extinction map from the previ- 
ously mentioned [Marshall et al.l (|2006l ) data. We utilise the 
UKIDSS GPS to ess entially improve the resolution of the 
[Marshall et al.l (|2006l ) data by analysing only stars that are 
situated directly along the same line of sight as the molecu- 
lar cloud. 

In Section [2j we describe the data used to derive the 
distance to each cloud, including the photometric data, and 
the data used to model the line of sight specific extinction- 
distance relationships. In Section [3] we discuss how we 
colour-colour and colour-magnitude select a specific sample 
of stars, of similar spectral type, to determine the line of 
sight extinction. In Section[4]we use the line of sight extinc- 
tion measurements to identify the spatial extent of molecular 
clouds. In Section [S] we derive, in detail, the distance to two 
example molecular clouds. We derive distances to a further 
19 molecular clouds and the results are presented in Sec- 
tion |B] We discuss the results and make comparisons with 
previous authors in Section [7] 



2 DATA 

2.1 The UKIDSS Galactic Plane Survey 



The U KIDSS Galactic Plane Survey (GPS, see lLucas et~ai] 
l|2008h ') covers, in the J(1.248Atm), H(1.631^m) and 



K(2.201/im) filters, the region of the Galactic Plane acces- 
sible by the United Kingdom Infrared Telescope (UKIRT); 
15° < ; < 107° and 141° < I < 230°, |b| < 5° and -2° < 1 
< 15°, |b! <2°. 

The GPS data are obtain ed from the WFCAM Science 
Archive (|Hamblv et al.l 120081 ) . The median 5a depths are 
J=19.77, 11=19.00 and K= 18.05 (Vega system) in the second 
data release of the GPS (| Warren et al.ll200'n ). The survey 
depth is spatially variable however, typically decreasing lon- 
gitudinally towards the Galactic centre and latitudinally to- 
wards the Galactic Plane. This occurs as the fie lds typically 
become more crowded towards these regions. iLucas et al.l 
(|2008l ) determine the modal depths in uncrowded fields to 
be J=19.4 to 19.65, H=18.5 to 18.75 and K=17.75 to 18.0. 

It is possible to use several different types of data qual- 
ity cuts, provided by the WFCAM Science Archive, to select 
only the most reliable data. To do this we remove sources 
with saturated pixels, remove blended objects using ellip- 
ticity cuts, astrometry cuts remove sources mismatched be- 
tween different filters, and finally photometric error cuts. 
The data cuts use d in this p aper arc almost identical to 
those presented in lLucas et al. (2008). The only alteration 
we make is a photometric error cut o f 0.03 mag in each 
photometric filter, whereas iLucas et al.l (120081 ) apply a 0.03 
mag colour cut. The use of such data quality cuts means 
that the photometric depths of the data used in this paper 
are reduced to J~18.0, H~16.8 and K~16.1. 



2.2 The Stellar Population Synthesis Model of 
the Galaxy 

We utilise the Stellar Populat ion Syn thesis Model of the 
Galaxy, developed in Besangon (jRobin et al.in2003i ). to model 
the distribution of giant stars, related to both the distance 
and extinction, along specific lines of sight. If an appropri- 
ate extinct ion- distance relationship is used to describe the 
extinction along the line of sight, then as will be shown, 
comparisons between synthetic and real data can be made 
to determine the distance to molecular clouds. 

A complete description of the Besangon model inputs 
can be found in lRobin et al.l (|2003l ). however we summarise 
it here for completeness. The model contains four popula- 
tions of stars, thin disc, thick disc, bulge and spheroid. Each 
of the four populations are described by a star formation 
rate history, an initial mass function, an age or age range, 
metallicity characteristics, kinematics, a set of evolutionary 
tracks, and includes a white dwarf population. The extinc- 
tion is modelled, in terms of visual magnitudes per kilopar- 
sec (mag kpc^^), by a diffuse thin disc. It is also possible to 
insert discrete clouds with a specified Av and distance. 



2.3 Distance-Av relationships 

[Marshall et all (|2006l ') use 2MASS data and the Besangon 
Galactic model to map the Galactic interstellar extinction 
distribution in three dimensions. This has been done to over 
64,000 lines of sight, in the inner Galaxy, each separated 
by 15' (hereafter referred to as the M06 distributions). In 
this paper we use the M06 data to assess the extinction- 
distance relationships along specific lines of sight. Using the 
M06 data, the distribution of synthetic giant population is 
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modelled by using a set extinction model to describe the 
thin disc (mag kpc~^) and the insertion of discrete clouds 
at specific distances. 

The M06 data produce reliable 3D extinction maps and 
so provide good line of sight specific extinct ion- distance re- 
lationships. In some cases it is possible to detect, and assign 
a distance to, large molecular clouds using the M06 data 
alone. We use the increased sensitivity of the UKIDSS GPS, 
in comparison to 2MASS, to improve the spatial resolution 
of the M06 data, identifying only stars that are situated di- 
rectly along the same line of sight as the target molecular 
cloud. 



2.4 The Galactic Ring Survey 

It is essential to know the position and spatial size of each 
molecular cloud in order to identify the st ars that are situ- 
ated d irectly along the same line of sight. iRathborne et al.l 
(|2009l ') have identified 829 mole cular clouds in the "C O J=l- 
Galactic Ring Survey fGRS I Jackson et al.ll2006l ). listing 
the Galactic position and the longitudinal and latitudinal 
FWHM of each cloud. The GRS covers a range of Galactic 
longitudes from 18° < 1 < 55.7° and |6| <1°. The GRS is 
fully sampled with a pixel size of 22". The processed data 
cubes have a Vlsr range from -5 to 135 km s~^ and a spec- 
tral resolution of 0.13 km s^'^ for 1 < 40°. For larger Galactic 
longitudes the Vlsr range covers -5 to 85 km s~^ and the 
spectral resolution is 0.21 km s^^. 

The GRS data serve two purposes in this paper. For 
comparison with the extinction maps produced in section 
14.31 and to identify the Galactic position and size of several 
clouds in order to determine their distances. 



3 ISOLATING RED GIANTS 

Unlike dwarfs stars, late type giant stars have lived long 
enough to have rotated around the Galaxy several times 
and have moved, in reference to the molecul ar clouds, from 
their original birth places in the spiral arms (jMarshall et al.l 
120061 ). For this reason it can be considered that late type 
giant stars are ubiquitous throughout the entire Milky 
Way. As they are infrared bright, numerous, have known 
colours and are found everywhere, they have been used 
by many authors to crea t e extinction maps of t he Galaxy 
(iLombardi fc Alvej l200ll : [Marshall et"ai] l2006l : iLombardil 
I2OO9I ). 

Using reliable photometry, it is possible to identify 
giant stars as they occupy known regions on colour- 
colour and colour-magnitude diagrams (CCDs/CMDs). 
Fig. [T] demonstrates this process with both Besangon 
and UKIDSS data. The synthetic Besangon data have 
been extracted from the same size and region of sky 
as the UKIDSS dat a, have been given a realistic ex- 
tinction distribution (jMarshaU et al.l 120061 ) and have been 
reddened along UKIDSS reddening tracks (jStead fc Hoard 
[2OO9I), and given realistic photometric errors. (The red- 
dening tracks used in this paper can be downloaded 
at www.ast.leeds.ac.uk/RMS/ReddeningTracks/). For these 
reasons the synthetic Besangon data can be considered a re- 
alistic representation of the real UKIDSS data. 

On a CCD, the most reddened stars (H-K>0.8) appear 



Table 1. CMD/CCD III extraction 



Data cut Besangon UKIDSS 



All 
CCD 

CMD 



6,591 
3,792 
1,339 



6,638 
3,329 
1,135 



to split into two separate branches of data. As the most red- 
dened stars have to be the most intrinsically bright, these 
two branches will be primarily composed of late giants and 
early dwarfs. As late giants are intrinsically redder than 
early dwarfs, the upper, and therefore redder, branch of data 
will be composed of giants, the lower will be composed of 

dwarfs. 

A GOIII stellar atmosphere model (jCastelli fc Kuruca 

[2OOJ) is progressively redden ed, using an extinction law of 
a=2.14 (jStead fc Hoarell2009f ). to create a reddening track, 
hereafter referred to as a GK04 GOIII reddening track, that 
splits the CCD, plotted in Figs. [T] (a and b), into two halves. 
The upper half contains stars redder, and therefore of a later 
spectral type, than GOIII. The CK04 GOIII reddening track 
has been shifted down the y-axis by ~0.04 to include stars 
that appear intrinsically bluer than a GOIII but whose la 
errors would allow them to be considered as GOIII or later 
candidates. 

As the synthetic Besangon data can be considered an 
accurate representation of the real UKIDSS data, and the 
class and spectral type of each synthetic star is known, we 
can estimate the make-up of our final UKIDSS sample. Late 
giants have very similar colours to late dwarfs of the same 
spectral type and therefore the colour-colour cut data will 
be contaminated with a large number of late dwarf stars. As 
these stars are faint, they can not be observed through large 
amounts of reddening. In colour-colour space the late dwarfs 
blend with the late giants in the field. In colour-magnitude 
space however, there is a clear distinction between giant and 
dwarf stars. Fig. [T](c) contains the previously selected sam- 
ple of GOIII and later candidates. Highlighted in black and 
red, in the synthetic data, are the dwarf and giant stars re- 
spectively, the remaining blue points consist of subgiants, 
bright giants and supergiants. The dwarf and giant stars 
form two clumps of data on the CMD. A black line has been 
placed, by eye, approximately equidistant from each clump 
separating the CMD into two halves. Although a small num- 
ber of giants are positioned in the first half of the CMD, it is 
safe to assume that the second half of the CMD is free from 
dwarf stars. We select this second half as our final sample 
of giants, and repeat this halving of the CMD process with 
the real data in Fig. [T] (d) . 

Should this method be applied to regions away from the 
Galactic Plane, where the line of sight extinction is lower and 
therefore the distinction between the giant and dwarf stars 
is less apparent, it may be necessary to statistically separate 
the two populations. However as there is a clear distinction 
between populations for all regions studied in this paper, 
such a process is not required. 

Table [1] contains the number of sources that remain af- 
ter each cut. The Besangon data do not model stars with an 
infrared excess. Such stars are abundant in star forming re- 
gions and are positioned to the right of the bulk of the data 
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Figure 1. a (Top left): A CCD of Besangon data modified as in lStead fc Hoard 1 120091 ) to replicate tiie UKIDSS data eentred on the 
RMS source G48. 9897-00. 2992 (see b). Stars tliat are redder thian a GOIII (red points) liave been separated from tlie remaining data 
using a CK04 GOIII reddening track (black line), b (Top right): Same as (a) except using UKIDSS data centred on the RMS source 
G48. 9897-00. 2992, covering an area of 60'x6'. c (Bottom left): A CMD of the remaining Besangon data (red points in (a)). Dwarf and 
giant stars are displayed as black and red points respectively and form two clumps of data. A black line has been placed approximately 
equidistant from each clump separating the CMD into two halves. Data to the right of the black line have been selected as the final 
sample. The remaining red points are subgiants, bright giants and supergiants. d (Bottom right): A CMD of the remaining UKIDSS data 
(red points in (b)). Like in (c), a black line has split the CMD into two halves, data to the right of the black line have been selected as 
the final sample. 



on a CCD, as such they are removed during the first data 
cut. In order to replicate the UKIDSS data, the Besangon 
data are clipped at the minimum and maximum J, H and 
K UKIDSS magnitudes. Of the final 1,339 synthetic stars, 
1,247, of a possible 1,626 from the total data set, are giants. 
Of the 1,247 giants in this final subset, 1,235 have a spectral 
type of GO or later. Therefore the reddest giants, spectral 
type GO or later, make up over 92% of the final UKIDSS 
sample. 

As shall be discussed in the following section, we dered- 
den the final UKIDSS subset to an intrinsic giant locus on a 
CCD, created using CK04 stellar models, in order to derive 
distances to the molecular clouds they surround. As over 
half of the 1,247 giants in the final Besangon subset have a 
spectral type between KO and Kl, we use a KOIII reddening 
track to deredden the UKIDSS data. Using a KOIII redden- 
ing track, created with the extinction law of lStead fc Hoard 



(|2009|), we calculate an approximate extinction-colour rela- 
tionship, Av~0.074(H-K). For each star a value of Aj is 
determined directly and converted to a valu e of Ay using 
the ra tio Aj/Av=0.2833, calculated using the lCardelli et al.l 
([1989|) extinction curve with Rv=3.1 (obtained from the Tri- 



legal website http://stev.oapd.inaf.it/cgi-bin/trilegaH ). It is 
noted in the literature that the value Rv=3.1 may not be 
appropriate for molecular clouds, however for the purpose 
of this paper it provides only a scaling between Ax and Ay 
and does not effect the errors of any derived results. 



4 MOLECULAR CLOUD IDENTIFICATION 

To illustrate the dereddening process and how the presence 
of molecular clouds along lines of sight can be detected, we 
obtain two samples of stars extracted from two different. 
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15'xl5' areas of the sky. The first area has been centred 
on the RMS source G048.9897-0.2992 and therefore has a 
molecular cloud along the line of sight. The second area is 
centred on the Galactic coordinates G048. 9897-0.7000. The 
area was visually selected from a region in the GRS that 
showed no obvious signs of ^'^CO, and therefore the line of 
sight extinction is likely to be free from dense molecular 
clouds. 



4.1 The dereddening method 

The stars along each line of sight have been colour-colour 
and colour-magnitude selected, as described in the previous 
section, to obtain a sample of stars primarily composed of 
red giants. We use a CK04 KOIII reddening track to dered- 
den each star in the sample created using the average ex- 
tinction power law of lStead fc Hoard (|2009l ). However to ac- 
count for the error created by dereddening stars of differing 
spectral types along a KOHI reddening track, instead of red- 
dening tracks specific to the spectral type of each star, we 
increase the error in the average extinction power law to 
Q=2.14lQgg. Each star is dereddened from its position on a 
CCD, along the KOHI reddening track, to the point of inter- 
section on the giant locus. As the KOHI reddening track has 
been created by convolving the progressively redder spectra 
of a CK04 KOIII stellar model through the UKIDSS filters, 
the length of the reddening track traversed directly relates 
to the amount extinction each star suffers. The photometric 
errors relate to the error in the point of intersection with 
the giant locus, and therefore relate to the error in the de- 
termined extinction. 



4.2 Distance-Av plots and Av histograms 

A distance-Av plot for the region G48.9897-0.2992, using 
the ~6,000 synthetic stars from Fig. [T] (a), is presented 
in Fig. [2] (a). It contains an extinction distribution model 
(EDM) created to follow the M06 data, extracted from the 
15'xl5' tile that is centred spatially closest to the RMS 
source. At ~5 kpc the EDM possesses a sharp rise in Ay. 
The position of this rise is consistent with the presence of a 
molecular cloud at the RMS kinematic distance of 5.2 kpc. 

Stars can also be binned into small groups representing 
small increments of extinction to construct Ay histograms. 
When plotting an Ay histogram of the synthetic data cen- 
tred on the region G48. 9897-0. 2992, the sharp rise in extinc- 
tion in Fig. [2] (a), at ~5 kpc, represents a fall in the star 
count over this particular range of Ay. In comparison, the 
same synthetic data, reddened using a smooth extinction 
distribution model of 1.4 mag kpc"'^, produces a steady rise 
and fall in star counts of the over a similar range in the 
Ay histogram (Fig. [2] (b)). If a distance- Ay plot could be 
constructed that contained precise distance and Ay mea- 
surements for each star in the field of view then, cissuming 
no stars are situated inside the molecular cloud, the molec- 
ular cloud would be represented by a vertical rise in the 
extinction distribution at the exact distance to the molecu- 
lar cloud. The corresponding Ay histogram would contain 
a sharp fall and rise in star counts, with a gap equivalent to 
the total Ay of the molecular cloud. 

Using only photometric data, it is possible to determine 



the amount of extinction a particular star suffers if accurate 
intrinsic colours are known. The amount of extinction suf- 
fered by a star is proportionate to how far the star moves, 
along a reddening track, on a colour-colour diagram (CCD), 
and can therefore be determined by dereddening individual 
stars along reddening tracks to their corresponding intrinsic 
colours. In this case the colour selected giants are dered- 
dened to a giant locus. The giant locus has been generated 
using UKIDSS filter profiles and the CK04 stellar models 
(|Stead fc Hoard I2OO9I ) . If we can accurately determine the 
amount of extinction suffered by every star along a partic- 
ular line of sight, then by creating Ay histograms the pres- 
ence of a molecular cloud can be detected by the drop in 
star counts. 



4.3 Extinction maps 

Fig- E] (a and b) contain Ay histograms of real 
stars extracted from the two previously mentioned re- 
gions G048. 9897-0. 2992, centred on an RMS source, and 
G048.9897-0.7000, centred on a "CO featureless region of 
space. As the Ay histogram that is centred on an RMS 
source will contain a molecular cloud along the line of sight, 
like the blue histogram in Fig. [2] (b) , there are two distinct 
peaks in the histogram, separated by a large dip in the source 
counts. The second region, that does not contain a molec- 
ular cloud, produces an Ay histogram with a steady rise 
and fall in star counts, and has a similar morphology to the 
black histogram in Fig. [2] (b). As noted, the presence of a 
molecular cloud in the former histogram creates two distinct 
peaks in the histogram. The far edge of the first peak there- 
fore represents the beginning of the molecular cloud. In the 
reverse manner, the near edge of the second peak therefore 
represents the end of the molecular cloud. To estimate the 
width of the gap, and therefore determine the Ay of the 
molecular cloud, two skewed Gaussians are fit to each peak. 
The respective la deviations are used to define the edge of 
each peak, and therefore the difference between them is de- 
fined as the Ay of the molecular cloud. In a similar fashion 
the far line of sight extinction is defined as the 3a deviation 
to the far side of the second peak. In Fig. [2] (b), as there 
is no molecular cloud and therefore no second peak, the 3a 
deviation is measured from a single skewed Gaussian. 

In this paper we have defined the far line of sight ex- 
tinction as the 3a deviation of all extinction measurements. 
This should be distinguished from the total line of sight ex- 
tinction, defined as the amount of extinction suffered along 
the line of sight towards the edge of the thin disk. Differences 
between the far and total line of sight extinction measure- 
ments can occur for two reasons. First, there may be a large 
amount of diffuse extinction behind the last detected cloud 
that would not fit well to a skewed Gaussian distribution. 
More importantly however, it could be that along the line of 
sight of the previously mentioned region G048. 9897-0. 7000, 
used to construct Fig.|3](b), there could be a molecular cloud 
with an optical depth large enough to prevent stars at the 
far side of the cloud being detected. Therefore the measured 
far line of sight extinction, that is dependent upon the depth 
of the photometry used, would be very different to the total 
line of sight extinction. However, the good agreement be- 
tween real data and synthetic data (see section 15.21 Fig. [Sjl 
suggests that the far and total line of sight extinction are 
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Figure 2. a (Loft); A distance-Av plot containing synthetic data showing the extinction distribution model (blue line) created to 
reproduce the M06 data (red error bars) centred spatially closest to the RMS source G048. 9897-0. 2992. The black line represents the 
smooth extinction distribution model of 1.4mag kpc~^. b (Right): A synthetic Ay histogram of each extinction distribution presented 
in the distance-Ay plot (Left panel). The blue data have a drop in star counts between Ay^l and ^^13, corresponding to the sudden 
increase in extinction beginning at ~4.5kpc in the distance-Ay plot. 



E 20 




200 




Figure 3. a (Left): An Av histogram of real stars extracted from a 15'xl5' region centred on the RMS source G048. 9897-0. 2992. We 
have defined the total line of sight extinction to the far side of the molecular cloud, as the la width (dashed blue line) on the left side 
of the skewed Gaussian (red line). (Right): An Ay histogram of a 15'xl5' region centred on an unremarkable region of space (see text) 
at G048. 9897-0. 7000. The histogram shows no presence of a molecular cloud. We define the far line of sight extinction (see text) as the 
3(T deviation of all extinction measurements (dashed green line) of the fitted skewed Gaussian (red line) . 



often the same. 

It is possible to repeat both of these processes in small 
increments for large fields to produce two different types of 
extinction map, far line of sight extinction (FLSX) maps and 
molecular cloud extinction (MCX) maps. The latter can be 
created by plotting Ay histograms and identifying which re- 
gions show the presence of molecular clouds. The Ay of each 
cloud can be measured and used to construct an extinction 
map, showing both the loca tion and size of an y clouds in the 
field of view. In comparison iLombardil ((20091), by measuring 
the colour excess of background stars, evaluate the extinc- 
tion at any location in the sky using a weighted average of 
extinction measurements to stars that are angularly close to 
a given location on the map. 

Using the above processes, both types of extinction 
maps have been created for a 2°x2° region centred at 



G49.5-I-0.0, containing the RMS source G048.9897-0.2992 
(The region was not centred closer to the RMS source due 
to incomplete UKIDSS data at 1<48.5). Using the deep 
UKIDSS data, we are able to produce Ay histograms for 
6'x6' tiles that are therefore over 6 times the resolution of 
the 15'xl5' M06 distributions. For each histogram we mea- 
sure the far line of sight extinction, and depending whether 
or not the histogram contains the presence of a molecular 
cloud, the Ay of the cloud is determined. Some histograms 
show the possible presence of more than one molecular cloud 
along the line of sight. In such cases we determine the far 
line of sight extinction to the far side of the most prominent 
molecular cloud. We analyse one 6'x6' region of sky and then 
microstep by shifting the field of view by 3', repeating the 
process for the entire 2°x2° region. The large region is then 
again split into 6'x6' tiles, however the extinction associated 
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Figure 4. a (Top): A FLSX map from the 2°x2° region centred at G49. 5+0.0. The 15, 18, 21 and 24 A^ levels are plotted in various 
shades of red. Blue contour linos represent the 11 Ay level. The RMS source is plotted as a red triangle. Blue + symbols represent 
regions where the extinction was too high to accurately determine the far line of sight extinction. The region of sky used to determine 
the M06 data is plotted as a dashed rectangle. A 9'x6' region surrounding the RMS source is plotted as a solid rectangle, b (Bottom): 
The same contour lines from (a) have been overlaid with a GRS integrated intensity image from the same region, created using the full 
velocity range available, from -5 to 85 km s~^. 
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Figure 5. a (Top): A MCX map from the 2°x2° region centred at G49. 5+0.0. The 4, 6, 8 and 10 Ay levels are plotted in various shades 
of red. The RMS source is plotted as a red triangle. Blue + symbols represent regions where the extinction was too high to accurately 
determine the Ay of detected clouds. The region of sky used to determine the M06 data is plotted as a dashed rectangle. A 9'x6' region 
surrounding the RMS source is plotted as a solid rectangle, b (Bottom): The same contour lines from (a) have been overlaid with a GRS 
integrated intensity image from the same region, created using a velocity range from 60 to 75 km s~ . 
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to each tile is taken from an average of the four surround- 
ing microstepped 6'x6' tiles. This process does not improve 
the resolution of the final extinction map, however it does 
smooth the extinction map, averaging out any anomalous 
Av measurements. 

Fig. U a contains a FLSX map from the 2°x2° region 
centred at G49. 5+0.0. Contour lines show the 11, 15, 18, 
21 and 24 Av levels. The RMS source G048.9897-0.2992 
(red triangle) is contained within the 21 Av level. Fig. 0] 
(b) contains the same extinction plot overlaid with a GRS 
integrated intensity image, from the same region, in order 
to draw a comparison between the positions of molecular 
clouds and the regions of high extinction traced out by the 
contours. The GRS image has been created using the full 
velocity range available from -5 to 85 km s^^. The contour 
lines have mapped out the bulk of the extinction shown in 
the GRS data, and from this we can estimate the extent of 
the molecular cloud that contains the RMS source. There 
are two other regions surrounded by Av=21 contour lines, 
at ~G48. 8+0.1 and ~G50.4-0.4, in both cases there do ap- 
pear to be other large groups of molecular clouds at these 
positions. 

Fig.[5]a contains a MCX map from the 2°x2° region cen- 
tred at G49. 5+0.0. Contour lines show the 4, 6, 8 and 10 Av 
levels of the measured depth of the molecular clouds. The 
RMS source G048.9897-0.7000 (red triangle) is contained 
within the 8 Av level. Fig. [S] (b) contains the same extinc- 
tion plot overlaid with the same GRS integrated intensity 
image, plotted in Fig. |4](b), however only the velocity range 
60 to 75 km s~^ has been considered. This range corresponds 
to the approximate 3a width of the measured kinematic ve- 
locity of the molecular cloud associated with the RMS source 
(|Urauhart et al.l 12008 '). Like in Fig. g] (b) the contour lines 
have mapped out the bulk of the extinction shown in the 
GRS data, and as before we can estimate the extent of the 
molecular cloud that contains the RMS source. Both types 
of extinction maps show very similar features, however in the 
MCX map, large areas are featureless due to the absence of 
molecular clouds. In all plots a 9'x6' region surrounding the 
RMS source has been overlaid. Stars within this region have 
been selected to further analyse in the following section. 



5 THE DISTANCE DETERMINATION OF 
MOLECULAR CLOUDS 

5.1 The distance-Av relationship 

Fig. [6]contains an Av histogram of real stars extracted from 
around the previously mentioned 9'x6' region surrounding 
the RMS source G048. 9897-0.2992. The morphology of the 
Av histogram is similar to the synthetic Av histogram, cre- 
ated using the EDM relating to the M06 data, presented in 
Fig. [2](b - blue histogram). As both histograms are similar 
it indicates that the M06 data are an accurate representa- 
tion of the extinction distribution along this particular line 
of sight. However there is an important feature worth not- 
ing, at Av~12 there is a drop in source counts creating a 
gap in the histogram, occurring due to the presence of the 
molecular cloud along the line of sight. If we estimate the 
range of Av along the histogram that the cloud covers, and 
we have an accurate line of sight specific extinction-distance 




Figure 6. An Av histogram of real stars from a 9'x6' region 
centred on the RMS source G048. 9897-0.2992. The near and far 
sides of the molecular clouds are marked with dashed blue lines, 
derived from the Icr deviations of the Gaussians (solid blue lines) 
fitted to the red bins. 



relationship, then it is possible to determine the distance 
to the cloud. The distance and the error in the distance to 
the cloud are derived by first measuring the Av values repre- 
senting the beginning and end of the molecular cloud. These 
Av measurements are then converted to distances using the 
extinction-distance relationships provided by the M06 data. 
The distance to the cloud is determined as the average of 
these two measurements, and the error, the difference be- 
tween them. The MOB error bars in Av are also considered 
and from this, we estimate the distance to the RMS source 
to be D = 5.5+0.4 kpc. This value is consistent with the 
RMS kinematic distance of 5.2+0.2 kpc. 

To calculate the certainty of detection of a molecular 
cloud in an Av histogram above the noise, a goodness of 
fit statistic is calculated as the square root of the summed 
squared residuals, presented in the following equation: 



^. Et^io^-E.r 



(1) 



where Oi is the actual star count and Ei is the expected star 
count, and n is the number of bins in the histogram. The 
goodness of fit statistic is calculated for both Gaussians, fit- 
ted to each peak in the histogram. A single Gaussian is then 
fitted to the same, combined range of data across both pre- 
vious peaks. A larger goodness of fit statistic is determined 
when fitting one Gaussian, thereby suggesting a poorer fit, 
than when using of two Gaussians. 



5.2 Measuring the Av gap 

In order to test the reliability of using two Gaussians to 
measure the Av of a molecular cloud on an Av histogram, 
we have measured the dip in source counts in 2500 different 
synthetic Av histograms. Each synthetic histogram has been 
created using a smooth EDM (EDMl) of 1.5 mag kpc^^ and 
contains a cloud with Av=1.5 at 5.2 kpc. To model the ef- 
fect of cloud patchiness, we also combine a smooth EDM 
(EDM2) of 1.5 mag kpc"'^ without a synthetic cloud along 
the line of sight. The ratio of EDM2 to EDM 1 is allowed 
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Figure 7. An Ay histogram of synthetic stars has been used to 
test the reliability of using two Gaussians to measure the Ay of 
a molecular cloud (see text). Each skewed Gaussian is shown as 
a red line and the corresponding la widths are shown as blue 
dashed lines. 



Figure 8. An Ay histogram of real stars (red histogram) taken 
from the same region of sky that was used to create the Marshall 
data (see text). The synthetic histogram, containing the same 
number of sources as the real histogram, is presented in black. 



to vary between 0% and 40%. Larger ratios than 40% make 
it difficult to detect the synthetic cloud above the noise in 
the histogram. Realistic errors are given to the synthetic 
data, and from these 2500 measurements, we compute the 
gap to be Ai/=1.6±0.5, suggesting that method used can 
reproduce the initial value of extinction used to create the 
synthetic cloud. An example synthetic histogram, produced 
using 20% of data from the cloudless EDM (EDM2), is pre- 
sented in Fig. [7] 

Including this additional error in Ay yields a distance 
to the RMS source G048. 9897-0. 2992 of D = 5.6±0.5 kpc. 
The additional error in the gap measurement of ±0.5 mag 
has had little effect on the resultant distance and error. 
This is because the extinction-distance relationship towards 
this particular region (at ~5 kpc), provided by the M06 
data, is visually estimated from Fig. [5] (a) to be ~5.5 mag 
kpc^^. Therefore large increases in the error in Ay will 
only translate into small errors in the derived distance. To- 
wards regions where the extinction-distance relationship is 
shallower, the additional error in the gap measurement will 
translate into larger errors in the derived distance. All sub- 
sequent distances derived in this paper include the error in 
the extinction-distance relationship, and the ±0.5 mag error 
in the gap measurement detailed above. 

Although Fig. [5] does contain an expected sharp Av 
gap at Av~12, the Kv histogram begins to show a reduc- 
tion in star counts at Av~8, continuing until Av~12. It is 
possible, due to the low number statistics, that this larger 
range is entirely caused by the molecular cloud and should 
be the measured dip. The distance assigned to this dip is 
D = 5.1±0.7 kpc. This value is still consistent with both 
the RMS kinematic distance and the distance previously de- 
rived. The error is smaller than might have been expected 
when measured from such a large range in Av, this is due to 
the steep slope on the EDM reducing the range of possible 
distances. However, we do not believe that the molecular 
cloud is responsible for this broader Av gap. There are sev- 
eral molecular clouds, identified in the Galactic Ring Survey, 
within the vicinity of the molecular cloud that contains the 



RMS source. The majority of these clouds are at a similar 
kinematic distance to G048. 9897-0. 2992 and so may be re- 
sponsible for the broadening of the Av gap. As the MOB data 
are not centred on the RMS source, the molecular cloud con- 
taining the RMS source is unlikely to be entirely responsible 
for the extinction measurements, along the line of sight, used 
to create the M06 data. 

The particular M06 data set used to model the ex- 
tinction along the line of sight extinction towards the 
RMS source G048. 9897-0.2992, is centred on the coordi- 
nates G049. 00-00. 25 and obtained from a 15'xl5' region. We 
have dereddened sources extracted from the same 15'xl5' re- 
gion as the MOB data and compared, in Fig. [S] a real Av 
histogram with a synthetic Av histogram. The synthetic 
histogram has been created by taking a random sample of 
Besangon data equal to the number of real stars. Both his- 
tograms show a similar morphology, however the real his- 
togram no longer has a second deep dip at Av~12. It is 
only when we extract stars from a smaller region centred 
exactly on the molecular cloud that this drop appears. This 
suggests that the Av dip in Fig. [S] is entirely due to the 
molecular cloud centred on the RMS source, and as such, 
the distance to the RMS source is, D = 5.B±0.5 kpc. 

The synthetic Av histogram in Fig. [5] was constructed 
with synthetic stars as distant as 15 kpc. As the real Av 
histogram also covers a similar range of Av as the synthetic 
Av histogram, it can be considered that for this particular 
line of sight, the far line of sight extinction and the total 
line of sight extinction, previously discussed in section 14.31 
are the same. 



5.3 RMS-GRS cloud associations 

As previously mentioned, Urquhart e t al (in prep. 



have identified the molecular clouds from iRathborne et al.l 
(120091) that are associated with 292 RMS sources. As 



IRathborne et al.l (|2009l ) have listed the Galactic position and 
the longitudinal and latitudinal FWHM of each cloud, this 
information can be used in the same manner as discussed 
to derive distances, independent to kinematic methods, to 
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RMS sources. The FWHM measurement of each cloud can 
be used to ensure we are observing stars directly in front 
and behind each cloud analysed. 

From the 292 RMS-GRS associations identified by 
Urqhart et al (in prep.), 173 have available J, H and K 
UKIDSS GPS data. To detect the presence of a molecular 
cloud in an Av histogram, enough foreground giant stars 
must be observed to notice the reduction in star counts per 
Ay bin that signals the presence of a molecular cloud. There- 
fore there will be a lower limit to the derivable distances to 
molecular clouds. If 3 star counts per Ay bin are needed to 
detect the presence of the foreground giants, in the synthetic 
histogram in Fig. [8] the first Ay bin containing this number 
of sources occurs when Av~3.5. Using the distance-Av plot 
in Fig. O created with the MOB data, converts this Ay value 
into a distance of ~2 kpc. As Fig. [2] contains a very typi- 
cal extinction distribution for the lines of sight considered 
in this paper, we therefore consider this a lower limit to 
the distances we can derive using this method and exclude 
clouds that have a near distance lower than 2 kpc. Further- 
more, we exclude clouds within 27° of the Galactic Centre to 
avoid the long thin Galactic bar, since this is not included in 
the Besangon model data. As the RMS and GRS kinematic 
distances have been derived using different Galactic rota- 
tion models, there are certain sources that do not possess 
an RMS kinematic ambiguity but are flagged as having so 
in the GRS catalogue. These clouds have also been excluded 
reducing our sample size to 74. It is not possible to assign 
an upper limit to the derivable distances, as the cut-off will 
be dependent upon the total extinction suffered and there- 
fore photometric depth. Clouds at higher Galactic latitudes 
may suffer less extinction than those at lower latitudes and 
it will therefore be possible to 'see further' along these lines 
of sight. We have assigned an initial upper limit of 9 kpc, 
and due to this upper limit, we have 34 clouds in our final 
sample that have a near distance larger than 2 kpc and a 
far distance smaller than 9 kpc. 

When using the GRS FWHM measurement to identify 
stars along the same line of sight as a molecular cloud, the 
area observed is very important. A large star count is needed 
to detect the presence of the cloud in an Ay histogram, how- 
ever, assuming each cloud to possess a spherical shape, stars 
along the line of sight of the edge of the cloud will suffer less 
extinction than those along the line of sight of the centre. 
Furthermore, the larger the area the greater the chance of 
contamination from overlapping clouds. To avoid excessive 
noise in the resultant Ay histogram we need to limit the 
size of the area observed, yet keep it large enough to detect 
enough stars. The size of each region observed is done on a 
cloud-by-cloud basis and we typically observe between a 0.6 
and 1.0 (7 width of the cloud. Finally, an EDM is created 
using the M06 data set that is spatially closest to the centre 
of the molecular cloud. 



5.4 Solving the kinematic distance ambiguity 

The previous region surrounding the RMS source 
G048. 9897-0. 2992 is particularly simple to analyse for 
two reasons. First, from inspection of the MOB data, there 
appears to be only one, very sharp rise in extinction. For 
this reason the distance to the cloud, using Fig. [2] (a), 
can be visually estimated to be 5.0±0.5 kpc. Second, the 



RMS source has a kinematic velocity located at a tangent 
point in the Galaxy and therefore has no near/far distance 
ambiguity. 

As a second demonstration we have extracted data 
along the line of sight of the molecular cloud GRSMC 
G041.04-0.66. The cloud has been associated with the RMS 
source G041. 0780-00. 6365 Urqhart et al (in prep.) and has 
two RMS kinematic distances of 5.8±1.1 and 7.0±1.1 kpc. 
Fig. [5] (a) contains a distance vs Ay plot of the extinction 
distribution model used to reproduce the MOB data assigned 
to the cloud. Unlike Fig. [2] (a). Fig. [9] (a) does not contain 
a sudden jump in extinction at the position of any of the 
kinematic distances. Therefore the resolution of the M06 
data is too low to distinguish between the near and far 
RMS kinematic distances. However, the M06 data still 
provide a reliable description of how the extinction relates 
to distance along the line of sight. 

Fig. |5] (b) contains an Ay histogram of the dereddened 
giants extracted from a la region of the molecular cloud. 
The histogram has two clear dips in star counts at Ay~8 
and Ay~ll that each separate two distinct peaks. We use 
the extinction distribution to plot the kinematic distances 
on the Ay histogram. The RMS and GRS kinematic 
distances are plotted as green and red arrows respectively. 
We first measure the dip that appears to be associated 
with the kinematic distances, and determine the distance 
to be D = B.9±0.5. Goodness of fit statistical tests have 
been applied to each Gaussian which have suggested that 
the dip is a real molecular cloud detection. This distance is 
consistent with the GRS kinematic distance D = B.38 kpc 
and the far RMS distance, D = 7.0 kpc. However as the 
near RMS distance is D = 5.8 kpc, the GRS distance is 
close to the average of the two RMS distances. Differences 
between the RMS and GRS kinematic distances occur due 
to the choice of the Galactic rotation model used. 

As previously mentioned, there is a second, equally 
well defined dip in the star count in Fig. |9] (b), signalling 
the presence of a second molecular cloud. This second dip 
appears at Ay~8 and is situated close to the near RMS 
kinematic distance. The distance to this second cloud has 
been determined to be D = 4.8±0.7 kpc. Although this 
value disagrees with the near RMS kinematic distance at 
the 1(7 level, there still remains some level of ambiguity to 
the derived distance to this molecular cloud. However, upon 
spatial examination of the entire GRS catalogue, it appears 
the region of sky previously studied also overlaps with a 
second molecular cloud G041.04-0.26, with a GRS distance 
D = 4.7 kpc. As the RMS source is spatially nearer the 
GRS cloud G041.04-0.66, the assigned distance to the RMS 
source is D = B.9±0.5 kpc. 

We have repeated the method on the remaining 33 
molecular clouds. All UKIDSS data towards each molecular 
cloud have been colour-colour and colour-magnitude cut 
to select a majority population of late giants. Goodness of 
fit tests have been performed on each and the results are 
detailed in the following section. 



6 RESULTS AND DISCUSSION 

The final 34 molecular clouds have each been sorted into one 
of three groups depending upon the features present in the 
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Figure 9. a (Left): A distance-Ay plot containing tlic extinction distribution model (black line) created to reproduce tlic M06 data (blue 
error bars) along the line of sight of the GRS molecular cloud G041.04-0.66. The GRS and two RMS kinematic distances arc plotted as 
red and green arrows respectively, b (Right): An Ay histogram of real stars from a small region centred on the GRS molecular cloud 
G041.04-0.66 contains a dip in star counts at Ay~ll. Two Gaussians have been fitted to each peak separated by the dip. The la widths 
have been used to determine the distance to the molecular cloud. 



Table 2. Glouds where sufficient data arc present to examine near and far distances. 



GRS Cloud 


Near 


Far 


GRS 


D 


error 




(kpc) 


(kpc) 


(kpc) 


(kpc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 



GRSMC G031.39+0.29 


6.35 


8.16 


6.55 


6.6 


0.5 


GRSMC G032.09+0.09 


6.13 


8.29 


7.07 


7.2 


1.1 


GRSMC G039.34-0.31 


4.55 


8.59 


4.55 


5.3 


0.7 


GRSMC G040.34-0.26 


5.34 


7.63 


5.43 


4.5 


0.8 


GRSMC G042.14-0.61 


4.99 


7.62 


5.12 


4.3 


0.5 


GRSMC G042.44-0.26 


4.87 


7.68 


4.90 


5.1 


0.7 


GRSMC G043.34-0.36 


4.46 


7.91 


4.65 


4.1 


0.7 


GRSMC G045.14+0.14 


4.53 


7.48 


4.53 


4.4 


0.6 


GRSMC G050.24-0.61 


3.22 


7.66 


3.20 


3.2 


0.2 


GRSMC G050.84+0.24 


3.43 


7.32 


3.53 


4.5 


0.4 


GRSMC G052.79-0.56 


4.78 


5.47 


5.78 


5.3 


0.6 


GRSMC G054.14-0.06 


3.76 


6.21 


6.32 


7.0 


1.0 


GRSMC G054.39-0.46 


2.76 


7.20 


6.82 


6.8 


0.8 


GRSMC G029.89-0.06 


6.14 


8.60 


6.78 


6.5 


0.7 


GRSMC G030.29-0.21 


6.76 


7.90 


7.32 


7.2 


0.7 


GRSMC G03 1.04+0. 29 


6.16 


8.44 


6.65 


6.3 


0.8 


GRSMC G033.04+0.04 


5.30 


8.96 


8.73 


8.1 


0.9 


GRSMC G041. 04-0.66 


5.78 


7.04 


6.38 


6.9 


0.5 


GRSMC G043.89-0.81 


3.92 


8.33 


3.85 


4.4 


0.6 


GRSMC G052.79+0.29 


4.55 


5.70 


6.10 


6.9 


1.1 



subsequent Ay histograms. For each of the final 34 clouds we 
have created both a distance- Ay plot, containing the MOB 
data and an extinction distribution model (EDM), and an 
Ay histogram. 

The first group contains 20 clouds and the correspond- 
ing Ay histogram of each has a strong presence of one or 
more molecular clouds. We derive the distances of these 
molecular clouds in an attempt to distinguish between the 
near and far kinematic distances. Results are presented in 
Table [2] and the table is split into two subgroups. The first 
subgroup contains regions where only one obvious molecu- 
lar cloud is present in the Ay histogram, such as the cloud 
containing G048.9897-0.2992 (Fig. [6]). The second subgroup 



contains regions where multiple molecular clouds are present 
in the Ay histogram, such as the cloud G041. 04-0. 664 (Fig. 
^jp). In cases where there are multiple molecular clouds, we 
use the cloud that appears to be associated with either the 
near or far distance, and address the presence of the ad- 
ditional clouds in the Appendix by looking for overlapping 
clouds in the GRS. The columns of Table [2] are as follows: 
(1) the molecular cloud name; (2) the near RMS distance; 
(3) the far RMS distance; (4) the GRS distance; (5) and (6) 
the distance and error derived in this paper. 

The second group contains 10 clouds and the corre- 
sponding Ay histograms do not contain an adequate number 
of sources with large enough values of Ay to exclude the far 
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Figure 10. a: A distance-Av plot containing the extinction distribution model (black line) created to reproduce the M06 data (blue 
error bars) along the line of sight of the GRS molecular cloud G042.04-0.01. The GRS and two RMS kinematic distances are plotted as 
red and green arrows respectively, b: An Ay histogram of real stars from a small region centred on the cloud G042.04-0.01. As data are 
insufficient to exclude the far kinematic distance, but there is also no obvious cloud at the near distance, the cloud has been placed in 
the first subset of Group 2 (see text), c: Same as (a) for the cloud G049.39-0.26. d: Same as (b) for the cloud G049.39-0.26 except the 
data are insufficient to exclude both the near and far kinematic distances. For this reason the cloud has been placed in the second subset 
of Group 2. e: Same as (a) for the cloud G043. 3900-0. 66. f: Same as (b) for the cloud G043.39-0.66, except no obvious clouds are present 
in the Av histogram and so the cloud has been placed in Group 3. 



kinematic distance. Results are presented in Tableland the 
table is split into two subgroups. In some cases it is still pos- 
sible to exclude the near distance due to the apparent lack of 
a molecular cloud at the near distance. The first and second 
subgroups contain regions where we have and have not been 



able to exclude the near kinematic distance respectively, Fig. 
IIUI (b) gives an example of an Av histogram where data are 
insufficient to exclude the far kinematic distance, but also 
contains no obvious cloud at the near distance, thereby solv- 
ing the near/far ambiguity through exclusion. Fig. [TD] (d) 
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Table 3. Clouds where data are insufRcient to exclude far distances. 



GRS Cloud 



(1) 



Near 
(kpc) 

(2) 



Far 

(kpc) 

(3) 



GRS 

(kpc) 

(4) 



Exclude Near? 



(5) 



GRSMC G030.44-0.26 
GRSMC G041.34+0.09 
GRSMC G042.04-0.01 
GRSMC G043.19-0.51 
GRSMC G044.29+0.04 
GRSMC G044.34-0.21 
GRSMC G045.49+0.04 



6.62 


8.04 


7.30 


Y 


4.12 


8.64 


8.55 


Y 


3.87 


8.74 


8.60 


Y 


4.15 


8.24 


8.27 


Y 


4.33 


7.83 


8.02 


Y 


5.33 


6.90 


6.80 


Y 


4.71 


7.22 


7.45 


Y 



GRSMC G028.59+0.04 


6.16 


8.76 


6.60 


N 


GRSMC G037.69+0.09 


6.18 


7.20 


6.70 


N 


GRSMC G049.39-0.26 


4.99 


6.04 


6.82 


N 



Table 4. Clouds where data are insufficient to detect the presence of a molecular cloud. 



GRS Cloud 


Near 


Far 


GRS 




(kpc) 


(kpc) 


(kpc) 


(1) 


(2) 


(3) 


(4) 


GRSMC G032. 99+0.59 


5.91 


8.36 


8.05 


GRSMC G030.79-0.06 


5.77 


8.84 


6.22 


GRSMC G043. 39-0.66 


4.32 


8.01 


8.18 


GRSMC G043.49-0.71 


3.32 


8.95 


2.83 



gives an example of an Av histogram where data are insuf- 
ficient to exclude both the near and far kinematic distances. 
The columns of Table [3] are the same Table [2] except as fol- 
lows; (5) a flag denoting if the near kinematic distance can 
be excluded. 

The third group contains 4 clouds and the correspond- 
ing Av histograms do not show the significant presence of a 
molecular cloud at either the near or far kinematic distance. 
Results are presented i n Table [H The columns o f Table [4] 
are the same Table [21 iRoman-Duval et al.l ((20091) provide 
the ^'^CO luminosity (K km s~^ pc'^) of each cloud in the 
GRS which can be used to estimate the Av of each cloud. 
This is done as follows; an average '^^CO luminosity of each 
cloud is determined by estimating the area of each from 
the Galactic longitudinal and latitudinal FWHM and the 
GRS distance. This value is then converted to a ^"^CO col- 
umn clump density, in cm~ '^, using the equation N(^'^CO) — 
8.75 X 10"W(^^CO), from lSimoneTal , (,20011). Finally this 
value is converted to Ay u sing the equation Ay = 4.24 x 
10""N("CO) + 1.67 from iGoodman et all (|2009l '). We do 
not expect this computed value of Av to agree particularly 
well with the value of Av, used to determine the distance 
to the cloud, measured in each histogram. This is because 
the computed value of Av has been averaged over the en- 
tire molecular cloud and the resolution of the GRS data is 
many times greater than the approximate, 6'x6' tiles used 
to construct each Av histogram. Furthermore, the spatial 
substructure of each cloud has not been considered. 

Three of the clouds in group 3 have the three lowest Av 
values, derived from the ^^CO luminosities, in our sample of 
34 clouds. This provides a likely explanation as to why the 
molecular clouds could not be detected above the noise in 
the Av histograms. For the remaining cloud, G30.79-0.06, 



the reverse is true. G30. 79-0.06 has one of the largest Av 
values, derived from the ^'^CO luminosities, out of the sam- 
ple of 34 clouds. Despite this however, it failed the goodness 
of fit test as the line of sight extinction is too large, from 
the combination of both the cloud G30.79-0.06 and a second 
molecular cloud, as shall be discussed in the Appendix, to 
detect an adequate amount of stars at the near side of the 
cloud. Fig. [To] (f ) contains an example Av histogram from 
group three. Distance-Av plots and Av histograms for the 
remaining regions, from all three groups, are available with 
the online material. 

To draw a comparison between our results and both the 
RMS and GRS kinematic distances. Fig. II II contains a plot 
of the GRS distances, both RMS distances and the distances 
derived to each of the 20 clouds in Group 1, subtracted from 
the GRS distance. The grey line crossing the origin there- 
fore represents the GRS distance subtracted from itself. The 
open blue circles represent the near RMS distance and the 
closed blue circles represent the far RMS distance. The red 
error bars represent the 20 distances derived in this paper 
from group 1. 

The near/far distance ambiguity has been resolved for 
27 clouds out of the sample of 34 molecular clouds. Distances 
have been derived to 20 clouds that are independent to kine- 
matic methods. Of the 20 distances derived, 19 agree the 
GRS distance to within the 2g level, and 18 of the 20 agree to 
within the Icr level, however this does not take into account 
the error on the GRS distance. One cloud, G050. 84-1-0.24 
with a GRS distance D = 3.5 kpc and a distance derived 
in this paper D = 4.5±0.4, agrees only within the 3o- level. 
There is a second overlapping cloud with a GRS distance of 
5.3 kpc, however there is no obvious sign of a second molec- 
ular cloud in the Av histogram, furthermore the distance 
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Figure 11. For the 20 molecular clouds in Group 1, the GRS distances, both RMS distances and the distances derived to each cloud 
in this paper have been subtracted from the GRS distance. The grey line crossing the origin therefore represents the GRS distance 
subtracted from itself. The open blue circles represent the near RMS distance and the closed blue circles represent the far RMS distance. 
The red error bars represent the distances derived in this paper. 



derived in this paper correctly resolves the near/far ambi- 
guity. 

In many cases the GRS distance agrees within errors 
with one of the RMS kinematic distances. However, as dis- 
cussed previously, due to the difference in choice of the 
Galactic rotation models used to derive both the RMS and 
GRS kinematic distances, there are some differences be- 
tween the two sets of kinematic distanc es. Error bars for the 
GRS d istances have not been plotted as lRoman-Duval et al.l 
((20091) do not include specific errors in each measurement, 
however they will be very similar to the size of the error 
bars on the RMS distances. Taking this into account, all 
distances measured in this paper agree well with the GRS 
resolved kinematic distances. 

There are often large differences between kinematic 
distances and those derived by independent means. For 
example, there are differences of ~2 kpc between kine- 
matic distances and distances derive d using trigono met- 
ric parallax in the Perseus spiral arm (|Xu et al.l 120061 ). Of 
the 18 ma s ers pr eviously mentioned that were studied by 
iReid et al.l ([20091), 3 lie within the same Galactic longi- 
tude range used in this paper. The differences between the 
kinematic distances and trigonometric parallaxes of these 3 
masers, presented bv lReid et al.l J2009), are in each case less 
than ~0.4 kpc. This value is therefore smaller than both the 
average error of the distances derived in this paper, and the 
average error of the RMS distances. 

One final limitation to this method is that the most 
reddened stars are predominantly M class giants, and as the 
most reddened stars, they are required to probe the most dis- 
tant molecular clouds. It has been noted previously that the 
MOIII stars should be dereddened along specific MOIII red- 
dening tracks. This is even more important for the most red- 
dened stars wher e the difference betw een reddening tracks 
is very apparent (jStead fc Hoarell2009l ). Although the error 



caused by using a KOIII reddening track has been accounted 
for, the error created is a systematic one, affecting all results 
in the same manner. It may be possible however to use the 
BesauQon data to estimate the position on a CCD where 
the dominant spectral type changes and adjust the choice 
of reddening track accordingly. The error caused by the in- 
appropriate choice of reddening track is ~0.03 mag in each 
photometric band, therefore similar to the error cut applied 
to the UKIDSS data in this paper. 



7 CONCLUSIONS 

Of the 34 molecular clouds studied in this paper, we have 
directly resolved the near/far kinematic distance ambiguity 
for 20 clouds by deriving a distance to each cloud. These 
distances have been derived using near-infrared photometry 
and an existing line of sight specific extinct ion- distance re- 
lationship. Therefore the distances presented in this paper 
are independent to kinematic methods typically used to de- 
termine distances to molecular clouds. These distances have 
been used to resolve the kinematic distance ambiguity asso- 
ciated with each cloud, an d from comparison to the work of 
iRoman-Duval et al.l (|2009l ) , all 20 clouds have been resolved 
correctly. 

For the remaining 14 clouds that we were unable to 
derive a distance to, 3 of the clouds were too optically shal- 
low to detect above the noise in subsequent Ay histograms. 
One cloud had several high extinction clouds along the same 
line of sight and so, due to a lack of foreground stars, we 
were unable to constrain a distance to the cloud. Finally, 
the method failed for 10 clouds as the extinction along the 
line of sight was too large to detect stars at, and beyond, 
the far distance. However, for 7 of these 10 clouds, we could 
not detect the presence of a molecular cloud at the near dis- 
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tance, thereby solving the near/far ambiguity through ex- 
cl usion. In aU 7 occur r ences our resuhs agreed with those 
of iRoman-Duval et ahl ()2009f ). In total, 27 clouds had their 
near/far ambiguity resolved. 

Although it was not possible to observe stars beyond the 
far RMS kinematic distance assigned to 10 of the 34 molec- 
ular clouds, it may still be possible to extend our sample to 
include clouds further than the original 9 kpc cut applied. 
The majority of derived distances are below 7.5 kpc, how- 
ever the method fails when the extinction along the line of 
sight is too large to detect stars at, and beyond, the far dis- 
tance. For this reason the method failed on some clouds, as 
close as 6.6 kpc, and succeeded on others as far away as 8.7 
kpc. Success of the method depends on the total extinction 
along the line of sight. Ahhough UKIDSS GPS data would 
allow the detection of stars suffering up to Av~25 using the 
error cuts detailed in this paper, from inspection of the Av 
histograms presented in this paper we can determine dis- 
tances to clouds that suffer below A\a~18 along the line of 
sight, including the visual extinction of the cloud itself. It 
may be possible to increase the depth by relaxing the error 
cut in the J band, the photometric band that suffers from 
the effect of extinction the most, without contaminating the 
resultant Kv histograms significantly. Furthermore, for re- 
gions where the extinction along the line of sight exceeds 18 
magnitudes of visual extinction, it may still be possible to 
confidently exclude the near distance. 

In this paper the GRS data provide a way to obtain the 
Galactic location and spatial extent of each molecular cloud 
studied. However the GRS only covers a Galactic longitude 
range from 1=18° and 1=55.7°, between |6| <1°. To extend 
the work in this paper to regions outside of this range, molec- 
ular cloud extinction maps (section 14. 3p could be used to 
identify molecular clouds and measure their spatial extent. 
This process would work best at larger Galactic longitudes 
where the average line of sight extinction is typically lower, 
making it easier to detect the presence of molecular clouds 
above the background extinction level. We will be able to 
extend and improve the method to clouds in the southern 
hemisphere when th e deeper VISTA data become available 
iMinniti et al.ll2009l '). 
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APPENDIX A: REGIONS WITH MULTIPLE 
MOLECULAR CLOUDS 

Distance-Av plots and Av histograms for all the remain- 
ing regions are available with the online material in Figures 
lAll to IA8I In all cases the distance-Av plots contain the 
extinction distribution model (black line) created to repro- 
duce the M06 data (blue error bars) along the line of sight 
of each GRS cloud. The GRS and two RMS kinematic dis- 
tances are plotted as red and green arrows respectively. The 
Av histograms contain the same red and green arrows cor- 
responding to the GRS and both RMS kinematic distances. 
In cases where a molecular cloud was present in the Ay his- 
togram, the two Gaussians used to estimate the extinction 
of, and from that the distance to, each cloud are shown as 
blue curves. The la width of each skewed Gaussian has been 
indicated with blue dashed lines. 



Al G029.89-0.06 

The molecular cloud G029. 89-0.06 has RMS distances of 
6.1/8.6 kpc and a GRS distance of 6.8 kpc. The distance 
derived in this paper is D = 6.5±0.7. Fig. lAll appears to 
contain a second molecular cloud at Ay~12. The derived 
distance to this second cloud is D = 8.3±0.7. This value 
agrees with the far RMS distance at the la level and so a 
certain level of ambiguity remains. 

The extracted region overlaps with a second molecular 
cloud G029.79-0.21. This second cloud, G029. 89-00. 26, has 
a high ^"^CO luminosity and a GRS distance of 8.45 kpc, 
therefore providing a likely explanation for the presence of 
the second cloud. 



A2 G030.29-0.21 

The Ay histogram, created from data extracted from around 
the molecular cloud, G030.29-0.21, appears to contain two 
Av dips of similar size (Fig. IA1|) . Measuring the size of the 
second dip, situated at Ay~7, reveals a distance of D = 
4.6±1.0 kpc. This distance is inconsistent with both RMS 
distances, 6.8/7.9 kpc, and the GRS distance, 7.3 kpc, as- 
signed to this cloud. 

The extracted region overlaps with the molecular cloud 
G030.14-0.06 that has been associated with the RMS 
source G030. 3830-00. 1099. This particular RMS source has 
a near/far ambiguity of 5.4/9.3 kpc. As this far distance is 
larger than our 9 kpc selection criterion, it explains why the 
RMS source was not present in our final sample. The derived 
GRS kinematic distance of the molecular cloud is D = 5.6 
kpc and so the presence of two molecular clouds in the Ay 
histogram has been accounted for. 



A3 G031.04+0.29 

The molecular cloud G031. 04-1-0. 29 has RMS distances of 
6.2/8.4 kpc and a GRS distance of 6.7 kpc. The distance 
derived in this paper is D = 6.3±0.8. Fig. IA2I appears to 
contain a second molecular cloud at Av~9. The derived 
distance to this second cloud is D = 7.5±0.6, therefore ex- 
cluding the far RMS distance. 

The extracted region overlaps with several molecular 
clouds G030.89-H00.14, G031.24-)-00.09 and G031. 29+00. 09 



that have GRS distances of 7.28, 7.25 and 7.25 kpc respec- 
tively. All three clouds have high ^^CO luminosities provid- 
ing a possible explanation for the presence of a second cloud. 



A4 G033.04+0.04 

The molecular cloud G033. 04-1-0.04 has RMS distances of 
5.3/9.0 kpc and a GRS distance of 8.7 kpc. The distance 
derived in this paper is D = 8.1±0.9. The second molecu- 
lar cloud present in the Av histogram, at Av~7 (Fig. I A3] ), 
has a derived distance D = 6.3±0.8. This excludes the near 
RMS distance at the la level. There are four other molecu- 
lar clouds nearby, all with GRS distances of 7.1 kpc, which 
could account for the presence of a second molecular cloud 
in the Av histogram. 

A5 G043.89-0.81 

The molecular cloud G043. 89-0.81 has RMS distances of 
3.9/8.3 kpc and a GRS distance of 3.9 kpc. The distance 
derived in this paper is D = 4.4±0.6. Fig. IA5I appears to 
contain a second molecular cloud at Av~9. The derived 
distance to this second cloud is D = 5.3±0.6, therefore ex- 
cluding the far RMS distance. 

The extracted region overlaps with a second molecular 
cloud G044.34-00.81 that has a GRS distance of 4.73 kpc, 
therefore providing a possible explanation for the presence 
of the second cloud. 



A6 G052.79+0.29 

The molecular cloud, G052.79-f 0.29, has RMS distances of 
4.6/5.7 kpc and a GRS distance of 6.1 kpc. The distance 
derived in this paper is D = 6.9±1.1 kpc. Fig. IA7I appears 
to contain a second molecular cloud at Av~12. The derived 
distance to this second cloud is D = 5.4±0.5, similar to the 
far RMS distance. There is however an overlap with a second 
cloud, G052. 64-1-0. 14, at a distance D = 5.1kpc providing a 
possible explanation for the presence of the second cloud. 



A7 G028.59+0.04 

The molecular cloud G028.59-f0.04 has RMS distances of 
6.2/8.8 kpc and a GRS distance of 6.6 kpc. The Av his- 
togram does not contain an adequate amount of data to 
detect a molecular cloud at the far RMS distance and so 
it has therefore been placed in group 2. There is however a 
very large dip in source counts in the Ay histogram between 
Ay~6 and ~14 (Fig. lAlj ). The distance assigned to this dip 
is D = 5.3±0.9. The error is much smaller than might first 
be expected for such a large range of Ay, however upon in- 
spection of the EDM there is a very large increase in Ay 
beginning at ~4.5 kpc. Although the distance does agree 
with the near RMS source at the la level, there are several 
other overlapping clouds covering a range of distances from 
~3.5 to ~7.5 kpc. 



A8 G045.49+0.04 

G045.49+0.04 has been placed in group two as it is not 
possible to exclude the far RMS distance (Fig. IA6|I . It has 



18 



RMS distances of 4.7/7.2 kpc and a GRS distance of 7.5 
kpc. There is however evidence of a molecular cloud in the 
Kv histogram at Av~6. It has an assigned distance of D = 
4.6±1.4 kpc and so we are also unable to exclude the near 
distance. For this reason the cloud has been placed in the 
second subgroup of group two. There is however a second 
molecular cloud along the line of sight. The second cloud, 
G045. 14+00. 14, has a GRS distance of 4.5 kpc. This cloud 
could therefore explain the presence of the molecular cloud 
in the Av histogram, and so it may actually be possible to 
exclude the near distance. 



A9 G030.79-0.06 

The molecular cloud G030.79-0.06 has RMS distances of 
5.8/8.8 kpc and a GRS distance of 6.2 kpc. It has failed 
the goodness of fit test and so has been placed in group 3. 
There is however a very strong presence of a molecular cloud 
in the Av histogram, not associated with either kinematic 
distance (Fig. IA2]| . The associated distance to this dip is D 
= 10.0±1.2 kpc. There is a second RMS source, G030.9281- 
00.0308, with an associated molecular cloud, G030.94-0.06, 
along the same line of sight. The RMS distances are 4.3/10.3 
kpc and the GRS distance corresponds to the near distance 
at 4.4 kpc. This would appear to be the only confrontation 
with the results of lRoman-Duval et al.l (|2009f ). There is how- 
ever another very CO bright molecular cloud, G031.24-0.01, 
with a GRS distance of 11.55 kpc that overlaps the same 
region of sky. 

As the GRS distance to the original cloud studied, 
G030.79-0.06, is 6.2 kpc, and the location of the dip, should 
it have been detected, would have been at the very near 
edge of the histogram fFig. IA2|I . it is apparent that for this 
line of sight we are unable to detect molecular clouds at dis- 
tances much smaller than the GRS distance D = 6.2 kpc. 
The method failed because the line of sight extinction is too 
large along this particular direction. The molecular cloud 
G030.94-0.06 has a GRS distance of 4.4 kpc, and so it is 
likely that we are unable to detect the presence of this cloud 
in the Ay histogram. As there does not appear to be a third 
Av gap in the histogram, it is likely that this prominent Av 
gap, at Av~ll really is associated with the molecular cloud 
G031. 24-0.01, at a GRS distance of 11.55 kpc. The M06 data 
set used is centred at G030.75-f 00.00 which is ~30' off-centre 
from this third molecular cloud G031.24-0.01. Therefore it is 
possible that the M06 data do not represent the line of sight 
extinction centred on this molecular cloud. For this reason 
the extinction measurement may not have been accurately 
converted to a distance measurement. 

This paper has been typeset from a T^iJX/ I^T^gX file prepared 
by the author. 
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